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[t]  We  report  on  subauroral  plasma  density  troughs  observed  by  Defense  Meteorological 
Satellite  Program  (DMSP)  satellites  during  the  major  magnetic  storm  of  6  April  2000.  The 
troughs  were  embedded  within  strongly  irregular,  >500-km  wide  subauroral  convection 
streams,  composed  of  two  distinctive  parts.  The  poleward  part  coincided  with  strong 
wave  structures,  energetic  (ring  current)  ion  precipitations,  and  enhanced  vertical  ion 
flows.  One  or  several  narrow  density  decreases  were  present  here.  The  equatorward  part 
was  marked  by  a  density  trough  coincident  with  a  local  maximum  of  the  electron 
temperature  ~7000-9000  K.  Here  the  ion  temperatures  and  downward  vertical  ion 
velocities  were  ~2000  K  and  <20  m/s,  respectively.  The  mean  convection  velocity  was 
typically  <500  m/s.  We  develop  analytical  approximations  for  the  rate  coefficients  of  the 
charge  exchange  reactions  based  on  recent  laboratory  experiments.  These  have  been  used 
in  local  modeling  of  the  equatorward  F  peak  density  depletions.  Topside  steady-state 
density  profiles  were  evaluated  assuming  diffusive  equilibrium  in  a  given  Te  profile.  A 
scenario  for  the  transition  of  an  initial  density-height  profile  to  a  flat  profile  is  described. 

Substantial  agreement  between  the  modeling  results  and  DMSP  observations  indicates  that 
the  vibrational  mechanism  contributes  significantly  to  the  formation  of  high-re-related 
density  troughs.  INDEX  Terms:  243 1  Ionosphere:  Ionosphere/magnetosphere  interactions  (2736); 
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1.  Introduction 

[2]  During  magnetic  substorms  and  storms,  deep  plasma 
density  troughs  form  in  the  evening  sector  at  subauroral 
latitudes  in  the  topside  ionosphere  [e.g.,  Spiro  et  al.,  1979; 
Evans  et  al.,  1983;  Anderson  et  al.,  1991,  1993].  These  are 
often  associated  with  enhanced  streams  of  westward  plasma 
convection,  described  by  the  names  polarization  jets  (PJ) 
[Galperin  et  al.,  1974]  or  subauroral  ion  drifts  (SAID) 
[Spiro  et  al.,  1979;  Anderson  et  al.,  1991,  1993].  Recently, 
Foster  and  Burke  [2002]  suggested  the  name  Subauroral 
Polarization  Streams  (SAPS)  to  include  plasma  flow  events 
that  have  both  broad  and  narrow  extents  in  latitude.  Stable 
Auroral  Red  (SAR)  arcs  are  also  observed  within  stormtime 
troughs.  Their  630.0-nm  emissions  are  caused  by  elevated 
electron  temperatures,  Te  <  3000  K  [e.g.,  Newton  et  al., 
1974;  Rees  and  Roble,  1975;  Foster  et  al,  1994],  These  are 
explained  as  consequences  of  heat  conduction  due  to  the 
ring  current-plasmasphere  interactions  [e.g.,  Kozyra  et  al., 
1987;  Jordanova  et  al.,  1999;  Liemohn  et  al.,  2000], 
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[3]  The  ion  population  of  the  topside  F  region  below 
1000  km  is  dominated  by  0+.  As  their  recombination  rate  is 
very  small,  ionization  balance  is  mainly  determined  by  rapid 
recombination  of  NO+  and  02,  stemming  from  the  ion- 
molecule  reactions 


(Rl) 

o+  +n2  ^no+  +n 

(R2) 

o+  +  o2  ^  o+  +  o 

The  rate  coefficients  kx  and  k2  vary  with  the  energies  of 
interacting  species.  Experiments  in  drift  tubes  specified  the 
dependences  of  kx  and  k2  on  the  ion  translational  energy 
e,-  cx  Vd  [e.g.,  McFarland  et  al,  1973;  Albritton  et  al, 
1977],  Here  Vd  =  cE0  x  B/B2  and  E0  is  the  electric  field 
(c  is  the  speed  of  light).  On  the  other  hand,  it  is  known  that 
k\  and  k2  also  increase  with  the  level  of  the  vibrational 
excitation  of  N2  and  02,  respectively  [Schmeltekopf  et  al, 
1967,  1968;  Hierl  et  al.,  1997],  In  thermal  equilibrium 
the  vibrational  population  is  defined  by  a  Boltzmann 
distribution  with  a  vibrational  temperature  Tv. 
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[4]  The  coincidental  formation  of  plasma  density  troughs 
and  SAPS  structures  has  been  explained  through  a  combi¬ 
nation  of  enhanced  recombination  stemming  mainly  from 
the  reactions  (Rl)  and  (R2)  and  plasma  outflows  due  to 
frictional  ion  heating  [Newton  and  Walker,  1975;  Schunk  et 
al.,  1975,  1976;  Banks  and  Yasuhara,  1978;  Anderson  etal., 
1991;  Moffet  et  al,  1998;  Pavlov,  1998;  Pavlov  et  al,  2000; 
Pavlov  and  Foster,  2001;  Vlasov  and  Kelley,  2003],  We 
define  the  trough  depth  as  8 nm  =  ( n0  -  »min)/«0,  where  «min 
( n0 )  is  the  minimum  (background)  electron  density.  Note 
that  hnm  —*  1  when  «mjn  — >  0. 

[5]  Newton  and  Walker  [1975]  and  Moffett  et  al  [1998] 
found  that  making  allowance  for  vibrational  excitation  of 
N2  at  Te  ~  3000-4000  K  reduced  the  electron  density  in 
the  F  region  by  a  factor  of  £2  (i.e.,  8 nm  £  0.5).  Using 
the  updated  rates,  Pavlov  and  Foster  [2001]  showed  that 
vibrational  excitations  at  Te  ~  3000  K  lead  to  8 nm  ~  0.65 
in  the  daytime  ionosphere  at  ~300-350  km.  Schunk  et 
al.  [1975,  1976]  showed  that  enhanced  recombination 
caused  by  intense  electric  fields  >0.1  V/m  yields  8 nm  a 
0.65-0.75  at  >300  km.  Note  that  E0  ~  0.1  V/m 
corresponds  to  horizontal  convection  velocities  VH  = 
Vd  ~  2  km/s.  In  addition,  Schunk  et  al.  [1975]  states 
that  the  presence  of  strong  vibrationally  excited  N2  with 
Tv  ~  4600  K  is  equivalent  to  the  effect  of  ~0.1  V/m 
electric  fields.  Overall,  these  investigations  agree  that  the 
vibrational  mechanism  cannot  solely  account  for  the 
trough  formation.  Anderson  et  al.  [1991]  suggested  that 
ion  outflows  with  mean  velocities  ~500-800  m/s  due  to 
frictional  ion  heating  are  the  main  cause  of  deep  troughs 
8«m  >  0.75-0.9. 

[6]  The  commonly  accepted  scenario  for  SAPS-related 
troughs  relies  on  plasma  streams  with  large,  typically  VH  > 
2  km/s,  velocities  and  strongly  elevated  ion  temperatures. 
Indeed,  usually  Te  enhancements  are  either  less  noticeable 
or  at  times  ~7).  For  example,  Moffet  et  al.  [1998]  presented 
observations  from  the  DE-2  satellite  at  altitudes  ~456,  415, 
and  315  km  of  three  troughs  with  VH  ~  2.3,  2,  and  4  km/s 
and  7]  (Te)  ~  4000  (3000),  4500  (4500),  and  5000  (2500)  K, 
respectively. 

[7]  However,  Forster  et  al  [1999]  and  Mishin  et  al. 
[2003]  reported  observations  of  unusual  subauroral  troughs 
during  the  21  March  1990  and  6  November  2001  magnetic 
storms.  The  first  one,  with  bnm  ~  0.67  at  ~700  km  altitude 
and  56.6°  magnetic  latitude,  was  encountered  in  the  dusk 
sector  (at  18.5  MLT)  during  a  nearby  pass  of  the  Active 
satellite  above  the  Millstone  Hill  radar.  The  electron  tem¬ 
perature  increased  to  7200  K  within  the  trough,  while  the 
ion  temperature  remained  <2000  K,  comparable  to  that  in 
the  adjacent  midlatitude  region.  Simultaneous  observations 
from  the  Millstone  Hill  radar  agreed  with  the  satellite 
observations  and  showed  that  the  density  depletion  and 
enhanced  Te  ~  3500  K  at  >400  km  altitude  were  embedded 
within  a  SAPS  with  VH  ~  0.7-1  km/s.  The  latter  agreed 
well  with  measurements  of  VH  ~  1 .2  km/s  at  ~830  km  by 
the  DMSP  F8  satellite  flying  ~50  min  earlier  and  later  than 
the  Active  pass.  Forster  et  al.  [1999]  also  noted  the 
presence  of  unusually  high,  ~1.5  x  102  cm-3,  molecular 
ion  densities  within  the  trough  at  ~700  km. 

[8]  The  6  November  2001  SAPS  event  was  observed  by 
the  DMSP  FI 5  satellite  at  0423  UT,  after  a  sharp  increase  of 
the  AE  index  from  ~500  to  ~2000  nT  following  a 


northward  turning  of  IMF  Bz  [ Mishin  et  al.,  2003],  Within 
a  6 nm  ~  0.75  trough  at  ~830  km,  Te  reached  ~9000  K, 
while  7]  ~  2000  K,  VH  ~  0.5  km/s,  and  a  downward  vertical 
flow  velocity  <20  m/s.  Significantly,  the  density  minimum 
occurred  ~1 00  km  equatorward  of  a  strong  SAPS  wave 
structure.  Here  the  term  “strong”  indicates  that  the  peak  to 
peak  variability  of  VH  oscillations  was  of  the  same  order 
or  even  exceeded  the  mean  convection  velocity  ( VH ). 
Enhanced  fluxes  of  precipitating  30  keV  ions,  likely  the 
low-energy  portion  of  the  ring  current  population,  were 
collocated  with  the  wave  structure  [Mishin  et  al.,  2003]. 
Apparently,  this  event  reflected  enhanced  ring  current- 
plasmasphere  interactions,  initiated  by  freshly  injected  par¬ 
ticles  [e.g.,  Anderson  et  al,  1993;  Gamer  et  al.,  2004], 
However,  a  special  study  is  needed  to  determine  how 
typical,  with  regards  to  the  trough  formation,  these  events 
are. 

[9]  The  objectives  of  this  paper  are  twofold.  First,  we 
report  on  the  trough  events  associated  with  the  strong 
SAPS  wave  structures  during  the  6  April  2000  magnetic 
storm  (section  2).  We  found  that  all  strong  SAPS  wave 
structures  were  accompanied  by  highly  structured  plasma 
density  troughs.  Typically,  one  of  the  density  minima 
coincides  with  the  wave  structure,  while  another  is  dis¬ 
placed  by  ~100  km  equatorward  and  coincides  with  the 
maximum  of  Te.  In  section  3  we  present  analytical 
approximations  for  the  rate  coefficients  of  the  reactions 
(Rl)  and  (R2j,  derived  from  the  latest  experimental 
results  [  Viggiano  and  Williams,  2001].  These  have  been 
used  in  local  modeling  of  the  equatorward  F  peak  density 
depletions,  performed  in  section  4.  To  compare  the 
modeling  results  with  the  DMSP  observations,  the  topside 
density  profiles  above  ~500  km  were  evaluated,  assum¬ 
ing  diffusive  equilibrium  in  a  given  Te  profile,  likewise 
Schunk  et  al.  [1975,  1976], 

2.  Observations 
2.1.  Instrumentation 

[10]  DMSP  satellites  are  three-axis  stabilized  spacecrafts 
that  fly  in  circular,  Sun-synchronous  polar  (inclination 
98.7°)  orbits  at  an  altitude  of  ~840  km.  The  geographic 
local  times  of  the  orbits  are  either  near  the  1800-0600 
(F13)  or  2100-0900  (F12,  14,  15)  meridians.  Owing  to 
the  offset  between  the  geographic  and  geomagnetic  poles 
DMSP  satellites  sample  wide  range  of  magnetic  local 
times  (MLT)  over  the  course  of  a  day.  The  ascending 
nodes  of  DMSP  orbits  are  on  the  duskside  of  the  Earth. 
Thus  the  satellites  move  toward  the  northwest  in  the 
evening  LT  sector.  Each  satellite  carries  a  suite  of  sensors 
to  measure  (1)  fluxes  of  precipitating  electrons  and  ions  in 
the  energy  range  between  30  eV  and  30  keV  (SSJ4);  (2) 
the  densities,  temperatures,  and  drift  motions  of  iono¬ 
spheric  ions  and  electrons  (SSIES);  and  (3)  perturbations 
of  the  Earth  magnetic  field  (SSM).  We  do  not  use  the  SSM 
data  in  this  study. 

[11]  SSJ4  sensors  consist  of  four  detectors,  one  high- 
energy  (1-30  keV)  detector  and  one  low-energy  (30- 
1000  eV)  detector  for  each  of  the  particle  types  [Hardy  et 
al.,  1984].  Nineteen  point  ion  and  electron  spectra  are 
returned  once  per  second.  SSIES  consists  of  (1)  an  ion 
drift  meter  to  measure  the  horizontal  ( VH )  and  vertical  (Vv) 
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Table  1.  The  6  April  2000  SAPS  Wave  Events 


Sat  Number 

Event 

UT 

MLat 

MLT 

FU 

11/23 

2330 

49° 

1930 

F12 

12/23 

2346 

-45° 

2042 

F14 

14/23 

2307 

-47° 

2042 

F15 

15/22 

2229 

49° 

2124 

F15 

15/23 

2344 

-46° 

2118 

cross-track  components  of  plasma  drift  within  the  range  of 
±3000  m/s  and  a  one-bit  resolution  of  12  m/s  for  ambient 
ion  densities  greater  than  5  x  103cm-3;  (2)  a  retarding 
potential  analyzer  to  measure  ion  temperatures  (7)),  com¬ 
position,  and  the  in-track  component  of  plasma  drift  (Fj|); 
(3)  an  ion  trap  to  measure  the  total  ion  density  («,);  and  (4)  a 
spherical  Langmuir  probe  to  measure  the  density  ( ne )  and 
temperature  (Te)  of  ambient  electrons  [Rich  and  Hairston , 
1994].  The  plasma  drift  and  density  measurements  are 
sampled  at  the  rates  of  6  and  24  Hz,  respectively,  while  Te 
and  Tt  were  sampled  at  a  rate  of  0.25  Hz. 

2.2.  SAPS  Wave  Events  on  6  April  2000 

[12]  The  main  phase  of  the  6  April  2000  magnetic  super¬ 
storm  began  at  ~1650  UT  and  reached  maximum  epoch  at 
~2400  UT  [e.g.,  Huang  and  Burke,  2004].  Throughout  the 
main  phase,  a  few  short-lived  magnetospheric  disturbances 
occurred  at  ~2100,  2200,  2300,  and  2320  UT.  After 
2200  UT,  several  strong  SAPS  wave  structure  events  were 
observed  in  the  evening  sector  by  four  DMSP  satellites. 
Table  1  lists  the  universal  times,  magnetic  latitude  (MLat), 
and  magnetic  local  times  (MLT)  of  the  observed  structures. 
In  general,  all  of  these  events  resemble  the  6  November 
2001  event  [Mishin  et  al,  2003]  (hereafter  designated  as 
Nov6).  We  consider  strong  wave  structures  as  the  bench¬ 
mark  for  enhanced  ring  current-plasmasphere  interactions, 
initiated  by  magnetic  activity. 

[13]  Figure  1  provides  an  example  typical  of  strong  SAPS 
wave  events.  It  was  crossed  by  the  DMSP  F14  satellite  at 
~2306:30  to  ~2308  UT  (hereafter  event  14/23)  during 
a  southern  hemisphere  pass  in  the  evening  sector  at 
~2042  MLT.  The  top  two  panels  show  directional  differential 
number  fluxes  of  downcoming  electrons  and  ions,  respec¬ 
tively.  Horizontal  (positive  sunward)  and  vertical  (positive 
upward)  components  of  the  convection  velocity  are  shown  in 
the  third  and  fourth  panels,  respectively.  The  fifth  panel 
shows  variations  of  the  electron  and  ion  temperatures 
sampled  with  4-s  time  resolution.  Here  dots  and  asterisks  or 
diamonds  and  triangles  designate  the  data  obtained  from  the 
Air  Force  Research  Laboratory  (AFRL)  [e.g.,  Rich  et  al., 
2003]  or  from  the  DMSP  SSIES  data  distribution  website 
(available  at  http://cindispace.utdallas.edu/DMSP/)  at  the 
University  of  Texas  at  Dallas  (UTD).  The  bottom  panel 
shows  1-s  average  values  of  the  ion  density. 

[14]  Attention  is  directed  to  the  following  four  features  of 
the  data.  (1)  Significant  fluxes  of  precipitating  1-30  keV 
ions  were  detected  coincident  with  the  strong  wave  struc¬ 
ture.  These  fluxes  probably  represent  the  low-energy  por¬ 
tion  of  the  ring  current  distribution  function  [Smith  and 
Hoffman,  1974],  (2)  Fluxes  of  ~30—  1 00  eV  ions  increased 
up  to  £3  x  10s  cm-2  s-1  ster-1  eV-1  within  the  wave 
structure  and  disappeared  at  ~2307:30  UT.  The  flux  of 
suprathermal  (<100  eV)  electrons  increased  slightly  to 


~3  x  104  cm-2  s"1  ster-1  eV-1  within  the  wave  structure. 
(3)  Electron  and  ion  temperatures  were  irregular  within 
the  region.  Apparently,  there  is  good  agreement  between 
the  AFRL  and  UTD  values  for  Te.  On  the  average,  Te 
increased  within  the  wave  structure  up  to  ~8000  K, 
descended  to  a  local  minimum  ~5200  K  at  ~2307:35  UT, 
increased  again  to  a  local  maximum  ~7500  K  at 
~2307:50  UT,  and  then  gradually  decreased.  Tt  values 
calculated  with  the  AFRL  and  UTD  algorithms  show  sig¬ 
nificant  difference  in  the  regions  of  energetic  ion  precipita¬ 
tion  within,  and  poleward  of,  the  wave  structure.  The 
difference  is  less  than  20%  after  ~2307:30  UT.  (4)  The 
plasma  density  had  two  distinct  minima,  ~3  x  104  cm-3, 
near  2306:30  and  2307:35  UT,  respectively.  The  former 
coincided  with  an  increase  in  the  energy  and  intensity  of 
precipitating  ions  simultaneous  with  enhanced  vertical  ion 
flows  ~±500  m/s.  The  second  trough  was  detected  after 
~2307:05  UT,  eel. 8°  equatorward  of  the  wave  structure. 
The  density  decreased  gradually  as  did  the  energy  of  supra¬ 
thermal  ions  and  Te,  regardless  of  the  varying  ion  tempera¬ 
ture.  At  2307:30  UT,  n,  dropped  by  a  factor  of  ~2  at  the 
equatorial  edge  of  the  ion  flux.  At  this  time,  the  vertical 
component  of  the  ion  flow  was  downward  ~20-40  m/s. 
After  reaching  the  second  minimum,  n,  slightly  increased 
and  then  decreased  again  to  ~4  104  cm-3  at  ~2307:45  UT. 
It  steadily  increased  afterward.  During  this  period  the 
ion  temperature,  mean  convection  velocity,  and  downward 
vertical  flow  velocity  remained  ~2000  K,  <400  m/s,  and 

10  m/s,  respectively. 

[15]  Similar  behaviors  characterized  all  of  the  equator- 
ward  trough  events  studied  during  the  April  2000  storm. 
However,  sometimes  within  the  wave  structure  multiple, 
latitudinally  narrow  density  decreases  were  present.  To 
study  effects  of  electron  and  ion  heating,  we  only  consider 
data  acquired  equatorward  of  the  wave  structures.  Other¬ 
wise,  the  contribution  of  precipitating  energetic  ions  and 
vertical  plasma  flows  must  be  accounted  for. 

2.3.  Plasma  Density  Variations 

[16]  Letjus  consider  how  relative  plasma  densities 
8n,/«o  =  8«i  =  (n0  —  «,)/n0  vary  with  the  mean  (0+)  ion 
energy,  e/L'i?,  in  the  center  of  mass  system  of  interacting 
species,  that  is  0+  and  N2  (or  02) 


_  3  m,mi,2  (Tj  |  Tffj  |  m,yj\ 

2  m,  +  mj.2  V/”i  mi,2  3k  )'  ' 

Here  n0  stands  for  the  plasma  density  at  the  equatorward 
edge  of  the  SAPS,  Vd  =  cE0  x  B/B2,  k  is  the  Boltzmann 
constant,  and  ()  designates  average  values  over  the 
Maxwell-Boltzmann  ion  and  neutral  distributions.  Indices 
1  and  2  stand  for  N2  and  02,  respectively.  As  usual,  the 
effective  ion  temperature  is  defined  as  7’[1,2)  =  |e/!_’2).  To 
simplify  calculations,  we  take  m,  ~  m2  ~  2 m,  and 
translational  temperatures  Tx  =  T2  =  Tn-  1000  K,  yielding 
=  £/2)„  =  £,_„  and  7/1*  =  =  Tin  [cf.  St. -Maurice  and 

Ton-,  1978].  _  _  _ 

[17]  By  definition,  5n,  increases  as  n,  decreases  (5n,  — >  6 nm 
at  nt  — >  nm in).  If  the  ion  translational  energy  controlled  the 


3  of  9 


A10301 


MISHIN  ET  AL.:  STORMTIME  SAPS-RELATED  TROUGHS 


A10301 


Figure  1.  Example  of  SAPS  wave  structure  observed  by  DMSP  F14  at  ~2307  UT  on  6  April  2000 
(event  14/23).  Top  panels  show  electron  and  ion  energy-time  spectrograms  of  downcoming  directional 
differential  number  fluxes  in  cm-2  s-1  ster-1  eV-1  (logarithmic  scale).  The  third  and  fourth  panels  show 
vertical  and  horizontal  components  of  the  convection  velocity,  respectively.  The  fifth  panel  shows 
variations  of  the  electron  (dots  and  open  diamonds)  and  ion  (asterisks  and  open  triangles)  temperatures, 
as  calculated  using  AFRL  and  UTD  algorithms.  The  bottom  panel  shows  1-s  averaged  ion  densities.  The 
vertical  dashed  line  indicates  the  equatorward  boundary  of  the  electron  plasma  sheet. 


plasma  density,  6 «,■  should  increase_with  e,_„  [cf.  Schunk  et 
al.,  1975,  1976],  Figure  2a  shows  6«,  variations  versus  e,_„ 
for  the  12/23,  Nov6,  14/23,  15/22,  and  15/23  events  (data 
gaps  were  present  during  the  event  11/23).  In  general,  5 nt 
appears  to  be  independent  of  e,_„  in  all  events  save  perhaps 
that  of  12/23.  The  latter  is  the  only  event  where 
the  convection  velocity  at  the  equatorward  side  rose  to 


~1000  m/s  ( E0  ~  50  mV/m),  yielding  e,_„  >  4500  K  (see 
note  in  section  3.1).  _ 

[is]  Figure  2b  shows  6n,  variations  plotted  as  functions 
of  Te.  In  general,  the  density  increased  (6 n,  — >  0)  at  the 
equatorward  side  of  the  trough  in  step  with  electron 
temperature  decreases.  A  linear  approximation  yields 
bm  ~  6«0  +  a  ■  10~4  •  Te  for  Te  >  3000  K,  where 
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Figure  2.  Variation  of  the  relative  plasma  density  (a)  with 
the  mean  ion  energy  e,-_„  and  (b)  with  electron  temperature 
in  five  events  studied  (see  text). 


a  ~  1  and  1.25  in  the  northern  (events  15/22  and  Nov  6) 
and  southern  (events  12/23,  14/23,  15/23)  hemispheres, 
respectively.  For  the  12/23  event,  this  tendency  is  violated 
in  the  vicinity  of  the  wave  structure,  consistent  with  a 
raise  in 

3.  Rate  Coefficients  of  Ion-Molecule  Reactions 

[19]  Variations  of  the  rate  coefficients  k\  and  k2  with  the 
energies  of  interacting  species  have  been  studied  in  numer¬ 
ous  laboratory  experiments.  Experiments  in  drift  tubes 
specified  the  dependences  of  kx  and  k2  on  the  ion  transla¬ 
tional  energy  e,  at  rotational/translational  Tx2  and  vibration¬ 
al  Ty'2'*  temperatures  of  molecular  gas  ~300  K  [e.g., 
McFarland  et  al.,  1973;  Albritton  et  al.,  1977],  Experiments 
in  a  high-temperature  flowing  afterglow  device  (HTFA) 
allow  studies  of  the  rate  coefficients  kx  and  k2  at  thermal 
equilibrium  7}  =  TX  2  =  'fy'7'  [e.g.,  Viggiano  and  Williams, 
2001].  Analytical  approximations  for  k{  and  k2  derived  from 
the  recent  HTFA  experiments  [Hierl  et  al.,  1997;  Viggiano 
and  Williams,  200 1  ]  are  presented  below.  Hereafter,  we  use 


the  normalized  temperature  Tj_  =  7^[K]/300,  where  S;  =  i,  e, 
n,  etc. 

3.1.  0+  +  N2 

[20]  Schmeltekopf  et  al.  [1967,  1968]  determined  specific 
rate  coefficients  kl  v  for  vibrational  states  of  N2  from  V-  0 
to  11  at  7}  =  T\  =  300  K,  that  is  A^iA^ohooic  =  cv.  While 
c |~1,  the  rate  coefficients  for  V  >  2  increase  dramatically; 
for  example,  c2  ~  38,  c2  ~  85,  C4  ~  220,  and  cs  ~  270. 
Comparing  drift  tube  and  HTFA  data  in  the  range  of 
temperatures  from  300  to  1800  K,  Hierl  et  al.  [1997] 
established  that  within  the  experimental  uncertainty 
(~20%)  the  dependence  on  e,-  for  vibrationally  excited 
nitrogen  is  the  same  as  for  the  ground  state.  Thus  specific 
rate  coefficients  can  be  represented  as  kx ^e,)  =  cv  ■  fc10(e,). 
Correspondingly,  k\  is 

[V2](0)+^C^N(K)Y  W 

v>  \  ) 

Here  [jV2]w  designates  the  population  of  the  V  vibrational 
state  of  N2.  An  analytical  approximation  for  ki  0  was  derived 
by  St.-Maurice  and  Torn  [1978]  from  the  Albritton  et  al. 
[1977]  drift  tube  data  (7)  =  pm, V2d) 

1012  •  *1  o  ~  1 .5  -  0.6t,-  +  0.09t?  at  1  <  Tj  <  5.7 

(5) 

~  2.7-  1.15t,  + 0.15i?  at  5.7  <  T|  <  20. 

[21]  Figure  3  shows  the  HTFA  [Hierl  etal.,  1997]  (circles) 
and  increased  by  20%  drift  tube  (squares)  [ Albritton 
et  al.,  1977]  data  points  along  with  1.2  ■  kl0  (asterisks). 
Scaling  compensates  for  small  differences  in  the  data  sets 
that  are  well  within  the  quoted  error.  The  rate  coefficient  kx  is 


£1  —  k\p 


T,  (K) 


Figure  3.  Rate  coefficients  for  the  reaction  (Rl).  Circles 
and  squares  show  the  HTFA  and  increased  by  20%  drift  tube 
data  points,  respectively;  asterisks  show  1.2  •  kX  Q.  Solid, 
dashed,  and  dotted  lines  show  kx  (equation  (4))  calculated 
with  k\p=  1.2  ■  kx  0  and  with  the  Boltzmann  distribution  of 
n\r>,  assuming  the  vibrational  temperature  Tv=Th  2000,  and 
3000  K,  respectively. 
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Figure  4.  The  rate  coefficient  for  charge  transfer  (reaction 
(R2)).  HTFA  and  (scaled)  drill  tube  data  points  are  shown 
by  circles  and  squares,  respectively.  Triangles,  dotted, 
dashed,  dot-dashed,  and  solid  lines  show  analytical 
approximations  10“u  •  k2$  and  k2  for  T„  =  1200  K  and 
Tv  =  300,  1000,  and  3000  K,  respectively. 


calculated  with  kxfi  — »  1.2  •  kl  0  and  with  the  Boltzmann 
distribution  of  [AI2]W-  The  calculation  results  are  shown  by 
solid,  dashed,  and  dotted  lines  for  vibrational  temperatures 
f-y  =  T„  2000,  and  3000  K,  respectively.  The  analytical 
expression  fits  the  experimental  data  extremely  well  [cf. 
Pavlov,  1998].  Note  that  the  ground  state  rate  coefficient  ki  fi 
increases  rapidly  only  for  large  7)-  >  2000  K,  consistent  with 
the  event  12/23  in  Figure  2a. 

3.2.  O*  +  02 

[22]  Hierl  et  al.  [1997]  also  derived  the  dependence  of 
k2  on  f)?'.  However,  Viggiano  and  Williams  [2001] 
showed  that  the  contribution  of  the  rotational  temperature 
T2  to  the  reaction  rate  is  also  significant.  Plotted  as 
functions  of  the  total  energy  z,ot  =  1.57)  +  T2,  the  HTFA 
[Hierl  et  al.,  1997]  and  drift  tube  [ McFarland  et  al., 
1973]  data  points  appear  identical  at  small  values  of 
e,0,  <  0.2  eV  (Figure  5  of  Viggiano  and  Williams 
[2001]).  Ascribing  the  difference  between  the  two  data 
sets  at  higher  energies  to  the  contribution  of  the  first 
vibrational  level  [Viggiano  and  Williams,  2001]  and 
assuming  a  Boltzmann  distribution  for  the  02  vibrational 
population  [O^^,  one  arrives  at 

k2  ~  10-' 1  ■*2,o{[02](0)+[02](1)(log10  (t4/)  -2)} 

,  (  2.2  _  (  34.04  \\  ^ 

2’°"l(^r' +605expr^JJ 


(t«  =  1  ,5t;  +  t„).  These  approximations  are  shown  in 
Figure  4  by  triangles  (10" 11  •  k2fi)  and  the  solid  line  (k2). 
Again,  the  analytical  expression  fits  the  experimental  data 
well.  Note  that  by  virtue  of  the  method  of  definition,  the 
approximation  in  equation  (6)  is  numerically  close  to  that  in 
Figure  3  [Hierl  et  al.,  1997]  at  7]  <  1800  K,  lagging 
significantly  behind  at  higher  temperatures.  We  must  now 


compare  expectations  for  energy  dependencies  (equations 
(4)  and  (6))  with  the  experimental  results  in  Figure  2. 


4.  Discussion 

[23]  Unlike  E0,  the  Te  enhancements  within  the  SAPS 
rapidly  decrease  at  altitudes  below  300  km  [e.g.,  Newton 
and  Walker,  1975;  Foster  et  al.,  1994],  Vibrational  excita¬ 
tion  due  to  electron  impact  depends  strongly  on  Te  (see 
Figure  5).  Furthermore,  the  number  of  “vibrationally 
active”  electrons  in  the  high-energy  tail  of  the  distribution 
decreases  when  the  ratio  nJ[N2\  falls  below  ~10-3  [. Mishin 
et  al.,  2000],  Thus  we  anticipate  little  contribution  from  the 
vibrational  mechanism  at  altitudes  below  ~300  km.  Con¬ 
sequently,  plasma  density  profiles  caused  by  the  vibrational 
mechanism  should  differ  from  those  caused  by  strong 
electric  fields.  In  particular,  we  expect  a  flat  density  altitude 
profile  in  the  trough,  similar  to  that  found  in  Figure  6c  of 
Foster  et  al.  [1994]. 

[24]  Schunk  et  al.  [1975]  stated  that  Tv  ~  4600  K 
produces  the  same  density  profile  as  the  electric  field 
~0.1  V/m  or,  equivalently,  Tin  ~  4300  K.  Indeed,  using 
the  Schmeltekopf  et  al.  [1968]  specific  rate  coefficients  that 
were  scaled  to  rates  at  300  K  (see  section  3.1),  one  gets  that 
at  Tv^  4600  K  the  charge  transfer  rate  ~  k]  fi(T,„  =  4300)  ~ 
10  •  kifi(Tin  =  2000).  Using  the  updated  rate  coefficients 
from  Figure  3,  one  obtains  k^{Tin  =  2000,  7>  =  3000)  ~ 
kifi(Tin  =  4300).  This  means  that  the  required  level  of 
vibrational  excitations  is  significantly  lowered.  This,  for 
example,  may  explain  differences  between  the  results  of 
Newton  and  Walker  [1975]  and  Pavlov  and  Foster  [2001]. 

[25]  Trough  depths  in  the  topside  ionosphere  are  some¬ 
times  significantly  smaller  than  those  near  the  F  peak  [e.g., 
Anderson  et  al.,  1991].  However,  they  are  sometimes 
comparable.  For  example,  one  can  evaluate  from  Figure  6c 
[Foster  et  al.,  1994]  that  6 nm  ~  0.85,  0.94,  and  0.78  at 
altitudes  800,  400  (F  peak),  and  300  km,  respectively.  We 
do  not  have  access  to  plasma  parameters  near  the  F 
peak  during  the  topside  ionosphere  events  studied  here. 
On  6  November  2001  at  ~0424  UT,  the  CHAMP  satellite 
[Riegber  et  al.,  2002]  (available  at  http://op.gfz-potsdam. 
de/champ/)  detected  Te  =  4800  K  and  8 nm  ~  0.7  at  ~450  km 
in  the  southern  hemisphere.  This  event  was  simultaneous 
with,  and  magnetically  conjugate  to,  the  strong  wave- 
structure  event  from  DMSP  F15  [Mishin  et  al.,  2003].  The 
density/temperature  values  in  the  topside  ionosphere  and 
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Figure  5.  Rate  coefficient  for  the  vibrational  excitation 
(V  =  0  — >  V  >  0)  of  N2  by  thermal  electrons. 
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Figure  6.  Variation  of  the  ki/ki,o-ratio  with  time  calculated 
at  Tj  =  2000  K  with  various  combinations  between  Te  « 
3200,  3500,  4000,  or  6000  K  and  ne  =  (1,  2,  or  3)  x  10s 
cm-3,  designated  in  label:  3/3500  stands  for  ne  =  3  x  105 
and  Te  =  3500  K,  et  cetera. 
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near  the  magnetically  conjugate  F  peak  are  close  to  those 
observed  during  the  21  March  1990  event  by  the  Active 
satellite  and  Millstone  Hill  radar  [Forster  et  al.,  1999]. 

[26]  To  evaluate  ionization  balance  at  altitudes  at  and 
below  400  km,  one  can  use  the  local  approximation  [e.g., 
Schunk  et  al.,  1975].  Further,  from  Figures  3  and  4  it  is  clear 
that  at  7’  >  2000  K  the  effects  of  vibrational  excitations  are 
more  significant  for  kx.  Since  [A/2]/[02]  £  10  in  the  topside 
ionosphere  [e.g.,  Hedin,  1991],  to  first  order  one  needs  only 
to  account  for  reaction  (Rl)  at  Tv  >  2000  K.  The  increase  in 
the  vibrational  population  of  N2  can  be  evaluated  from 

Q 

-Q{t\v  tzPvQyVo  -kyt]y+L0r\y+l.  (7) 

Here  Qy  and  Lv  are  the  production  and  loss  rates  of 
vibrational  levels  r\y  =  [A2](KV[A2],  respectively.  The  total 
rate  Q  of  vibrational  excitation  by  (Maxwellian)  thermal 
electrons  is  calculated  using  the  cross  sections  from  Majeed 
and  Strickland  [1997].  Figure  5  shows  the  dependence  of  Ql 
ne  on  the  electron  temperature.  The  branching  ratios  for  V= 
1,  2,  3,  and  4  at  Te>  0.3  eV  can  be  estimated  as  (3K  ~  0.6, 
0.3,  0.1,  and  0.05,  respectively.  The  contributions  of 
suprathermal  electrons  via  electronic  triplet  states  and  of 
transport  processes  are  neglected.  At  altitudes  >300  km  the 
main  quencher  of  vibrational  states  is  atomic  oxygen  with 
the  rate  coefficient  [McNeal  et  al,  1974]  L0  ~  10“ 10  • 
exp(— 70/7,y3)  •  [O]  s-1  ([O]  represents  the  atomic  oxygen 
density).  Accounting  for  charge  transfer  yields  Lv  =  L0  + 
Cy  *  ^i,o  *  ti,  s 

[27]  The  vibrational  distribution  is  far  from  equilibrium, 
as  V  —  F exchange  ~5  10-13  •  [V2]  s-1  is  insignificant  at 
altitudes  >300  km.  To  evaluate  the  charge  exchange  rate  k\, 
the  population  of  the  first  five  states  (^oftr  =  1)  was  found, 
assuming  r\v>5  =  0.  Figure  6  shows  kx/kxfi  versus  time, 
calculated  at  7}  =  2000  K  with  various  combinations 
between  Te  «  3200,  3500,  4000,  or  6000  K  and  ne  =  (1,  2, 
or  3)  x  105  cm-3. 

[28]  After  ~5-10  min  the  production  of  NO+  due  to 
reaction  (Rl)  intensifies.  The  recombination  rate  of  NO+ 
ions  is  5  1/min  for  ne>  10s  cm-3.  Thus  the  density  decay 
rate  is  of  order  ~  kx  ■  [/V2]  S;  0.1  min-1.  As  the 


excitation  rate  Qv  oc  ne,  the  density  drop  affects  the 
variation  of  r[y  and  hence  must  be  accounted  for.  To 
first  order,  the  ionization  balance  can  be  described  by 
reaction  (Rl)  and 

^tne=k\  (0) «, (0)  [A2]  -  ar  ( ne  -  n,)  ne,  (8) 

assuming  ne=nt  +  [M?+]  and  «,{0)  =  ne( 0);  the  recombination 
rate  coefficient  ar  ~  4  x  10-7t7°'5  cm-3  s-1. 

[29]  Figure  7  shows  the  plasma  density  decrease  in  vibra- 
tionally  excited  gas.  This  was  calculated  with  Tin  =  2000 
(dotted  lines),  2500  (dashed  lines),  3000  (dash-dotted  lines), 
and  3500  K  (solid  lines)  for  Te  =  3500  and  4000  K  and  for 
the  MSIS  parameters  of  neutral  gas  at  ~300  and  400  km, 
respectively.  Apparently,  within  the  decay  time  the 
plasma  density  decreases  by  a  factor  of  ~6  to  ~10,  that 
is  &«,  ~  0.8 -0.9. 

[30]  Figure  8  (top)  shows  the  variation  of  the  plasma 
density  at  ~300  km  with  time  for  Tm  =  2000  K  and  Tc  = 
2000-5000  K.  The  relative  plasma  density  at  the  end  of  run 
versus  Te  is  plotted  in  the  bottom  panel.  The  depth  of  the 
density  depletion  clearly  increases  with  Te.  Note  that  the 
decay  time  is  shorter,  and  the  depletion  is  deeper  at  larger 
Tin,  consistent  with  the  results  of  Schunk  et  al.  [1975]. 

[31]  The  decay  time  increases  and  the  local  approxima¬ 
tion  breaks  in  the  topside  ionosphere  above  ~500  km.  To 
compare  the  modeling  results  with  the  DMSP  observations 
at  ~840  km,  magnetic  field-aligned  transport  processes 
(ambipolar  diffusion)  must  be  included.  However,  steady- 
state  topside  density  profiles  can  be  evaluated  assuming 
diffusive  equilibrium  in  a  given  temperature  profile  [cf. 
Schunk  et  al.,  1975].  The  0+  scale  height  is  H{0+ )  ~  K(Te  + 
Td/(mig\\),  where  g|j  is  the  field-aligned  component  of 
gravitational  acceleration.  Schunk  et  al.  [1975]  obtained 
H{0+ )  ~  300-350  km  due  to  the  large  7}  oc  Fq.  In  the 


0  10  20  30  40  50  60 

time  (min) 


Figure  7.  Density  depletion  at  ~300  (bottom)  and  400  km 
(top):  Tin  =  2000  (dotted  lines),  2500  (dashed  lines),  3000 
(dash-dotted  lines),  and  3500  K  (solid  lines). 
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Figure  8.  (top)  Variation  of  the  plasma  density  with  time 
at  300  km  for  Te  =  2000, . . .  5000  K;  Tin  =  2000  K.  (bottom) 
Variation  of  the  relative  plasma  density  at  the  end  of  run 
with  Te. 


events  under  consideration,  even  larger  values  of  H(0+)  ~ 
400-500  km  are  achieved  due  to  elevated  Te. 

[32]  Thus  nearly  constant  density  versus  height  profiles 
should  be  established,  yielding  large  6 nm  from  the  bottom  to 
the  top  of  the  trough  [cf.  Foster  et  al,  1994],  Given  that  in 
the  events  under  study,  Te  at  ~800  km  is  about  twice  as 
large  as  that  at  ~300-400  km,  one  can  compare  the  topside 
observations  (Figure  8)  to  the  F  peak  modeling  (Figure  2b). 
Good  apparent  agreement  between  them  indicates  that  the 
vibrational  mechanism  is  capable  of  forming  deep  density 
troughs  or  at  the  very  least  is  a  major  contributor. 

[33]  One  can  outline  the  dynamics  of  the  depletion 

process,  as  follows  (no  horizontal  transport  is  involved). 
As  soon  as  the  F  peak  density  decreases,  the  initial  diffusive 
equilibrium  breaks  down.  Plasma  in  the  vicinity  of  the  F 
peak  flows  into  the  density  hole  to  restore  the  force  balance 
along  the  flux  tube.  The  plasma  response  resembles  a 
rarefaction  wave  propagating  from  the  F  peak  and  can  be 
described  in  terms  of  field-aligned  ambipolar  diffusion. 
The  propagation  time  in  the  topside  ionosphere  is  of  the 
order  of  several  ts  ~  H(0+)/vs,  where  vs  ~  \/^(  Te  +  Tj)/m,. 
For  the  above  parameters,  one  gets  ts  ~  5-10  min.  As 
the  F peak  density  decay  rate  0.1  min-1  increases  with 

ne,  “fresh”  plasma  increases  the  “processing”  rate  upon 
arrival. 

[34]  The  depletion  of  a  whole  magnetic  tube  ends  when 
the  rates  of  plasma  field-aligned  transport  in  the  topside 
ionosphere  and  decay  in  the  depleted  F  peak  match.  This  is 
achieved  when  the  plasma  density  drops  to  several  104  cm-3 
and  the  density  gradient  scale  height  becomes  ~500  km. 
The  time  during  which  trough  plasma  achieves  the  steady 
state  is  of  the  order  of  several  tens  of  minutes.  This  scenario. 


as  depicted  in  Figure  9,  predicts  downward  topside  plasma 
flows  within  the  trough,  consistent  with  the  observations. 

[35]  Two  final  remarks  are  in  order.  First,  deexcitation  of 
excited  states,  which  might  become  important  at  high  Tv. 
was  not  included  in  equation  (7).  This  should  slow  down  the 
process  at  low  Te  <  3000  K,  shifting  the  beginning  of  the 
saturation  in  Figure  8  (bottom)  to  somewhat  higher  temper¬ 
atures.  Second,  for  a  given  electron  heating  source,  lower 
electron  densities  lead  to  higher  electron  temperatures.  Thus 
the  thermal  energy  balance  equation  must  be  included.  Such 
modeling  is  well  beyond  the  scope  of  this  observation-based 
study. 

5.  Conclusion 

[36]  We  report  on  subauroral  plasma  density  troughs 
associated  with  the  strong  SAPS  wave  structure  events 
observed  by  DMSP  satellites  during  the  major  magnetic 
storm  of  6  April  2000.  The  troughs  were  embedded  within 
strongly  irregular  >500-km  wide  subauroral  convection 
streams,  composed  of  two  distinctive  parts.  The  poleward 
part  coincided  with  the  strong  wave  structure,  energetic 
(ring  current)  ion  precipitations,  and  enhanced  vertical  ion 
flows.  One  or  several  narrow  density  depletions  were 
present  here.  The  equatorward  region  was  marked  by  a 
density  trough  coincident  with  a  local  maximum  of  the 
electron  temperature  ~7000-9000  K.  Here  the  ion  temper¬ 
atures  and  downward  vertical  ion  velocities  were  ~2000  K 
and  £20  m/s,  respectively.  The  mean  horizontal  convection 
velocity  was  typically  <500  m/s. 

[37]  Analytical  approximations  for  the  rate  coefficients  of 
charge  exchange  reactions,  derived  from  the  latest  labora- 


Vibrational  excitation 


lon/neutral  chemistry 
F  peak  depletion 

|  Ambipolar  diffusion 


Density  profile  flattening 


Figure  9.  Block-diagram  showing  a  scenario  of  the 
formation  of  the  high-7],  related  density  trough  (PS  stands 
for  plasma  sheet). 
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tory  experiments,  were  used  for  local  modeling  of  the  F 
peak  density  depletions  in  the  equatorward  regions.  The 
topside  density  profiles  were  evaluated  assuming  diffusive 
equilibrium  in  a  given  Te  profile.  Substantial  apparent 
agreement  with  the  observations  indicates  that  the  vibra¬ 
tional  mechanism  significantly  contributes  to  the  formation 
of  the  deep  density  troughs. 
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